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Abstract. The excess above 1 GeV in the energy spectrum of the diffuse Galactic gamma radiation, 
measured with the EGRET experiment, can be interpreted as the annihilation of Dark Matter (DM) 
particles. The DM is distributed in a halo around the Milky Way. Considering the directionality 
of the gamma ray flux it is possible to determine the halo profile. The DM within the halo has a 
smooth and a clumpy component. These components can have different profiles as suggested by N- 
body simulations and the data is indeed compatible with a NFW profile for the diffuse component 
and a cored profile for the clumpy component. These DM clumps can be partly destroyed by tidal 
forces from interactions with stars and the gravitational potential of the Galactic disc. This effect 
mainly decreases the annihilation signal from the Galactic centre (GC). In this paper constraints 
on the different profiles and the survival probability of the clumps are discussed. 
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1 Introduction 



From WMAP measurements of the temperature aniso- 
tropics in the Cosmic Microwave Background (CMB) 
in combination with data on the Hubble expansion and 
the density fluctuation in the Universe pQ we gather 
that Cold Dark Matter (CDM) makes up 23% of the 
energy of the Universe. The nature of the Dark Mat- 
ter (DM) is unknown, but one of the most promising 
particles is the "weakly interacting massive particle" 
(WIMP). Assuming that WIMPS are Major ana par- 
ticles they can annihilate each other and produce a 
large amount of secondary particles. For the determi- 
nation of the density distribution, the so-called DM 
halo profile of the WIMP particles, the gamma radia- 
tion is most important because it is not influenced by 
the magnetic field of the galaxy and points back di- 
rectly to its source. 

The observation of the diffuse Galactic gamma radia- 
tion of the Milky Way with EGRET showed an excess 
above 1 GeV in the photon energy spectrum. This ex- 
cess is different for various sky directions and can be 
interpreted as a WIMP annihilation signal [5] . There- 
fore it can be used to determine the DM halo profile. 
In section [5J we will describe the mechanism of the de- 
termination of the DM halo profile from the EGRET 
excess and explain how to calculate the DM annihi- 
lation flux of gamma rays. Then, after differentiating 
between diffuse and clumpy DM, we will introduce a 
survival probability for DM clumps as well as a ringlike 
substructure of DM within the Galactic disc. 
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2 Determination of the halo profile 

This directionality of the gamma ray excess depends 
on the shape of the DM halo. Therefore a fit to the 
signal in all different sky directions can be used to de- 
termine a possible halo profile. 

We divided the sky into four different latitude ranges 
(absolute values of latitude: 0-5°, 5-10°, 10-20°and 20- 
90°) and 45 longitude bins. Consequently, we obtained 
180 different sky regions. Strong point sources in the 
galaxy like CRAB and VELA were excluded because 
of their large contribution to the flux of the photons 
with a peculiar spectrum. The local annhilation rate 
is proportional to the square of the number density of 
DM particles. Therefore the gamma ray flux is calcu- 
lated with the line-of-sight integral 

Here 6/ is the branching ratio into the tree- level an- 
nihilation final state while dNf/dE is the differential 
photon yield for the final state f. The value of the 
WIMP mass M x is an essential parameter for the cal- 
culation of the gamma flux and can be obtained by 
fitting the shape of the photon energy spectrum in 
different sky regions. A WIMP mass in the range of 50 
- 70 GeV yields a good fit (see Ref. Furthermore, 
p x is the DM density profile of the halo and B\ is a 
local enhancement factor, the so-called "boost factor" 
which characterizes the increased WIMP annihilation 
rate in subhalos (DM clumps) in the DM halo of the 
Milky Way. These DM clumps are predicted by ana- 
lytical calculations [4] and numerical simulations [5] 
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from the inflationary-produced adiabatic density fluc- 
tuations. In Ref. [5] Eq. [Tj was used to determine the 
halo profile. However, for a clumpy DM distribution 
one expects the flux to be proportional to p x , c iump in- 
stead of p x c i ump because the clumps have a certain 
luminosity p c i,p- Therefore the flux is proportional to 
pd,p ■ p x ,cium P instead of p XjC i ump - The boost factor 
can in this case be written as [5] 

o , , I Px,cium P (r) dr g P (r) p c i,p{ r i n P ) fn s 

Hi = H . . r . (2) 

In Eq. [2] the parameter £p (r) characterizes the mass 
fraction of clumps in the DM halo, n p is the power 
index of primordial density perturbations. 
The first term in Eq. [2] represents the boost factor for 
the diffuse DM component while the second term is the 
boost factor for the clumpy DM. Assuming that the 
mass fraction in clumps is larger than a few percent the 
whole annihilation rate is dominated by the clumpy 
DM flux which can be approximated as 

$DM, clump 0^ 

(3) 

However, the expression fp(r, n p ) ■ /O c ;,p(r, n p ) is equal 
to a boost factor B, which is constant for all directions, 
multiplied with a factor P(r), which is a "survival 
probability" , averaged over all clump masses as calcu- 
lated in Ref. [6] . There it was assumed that clumps are 
completely disrupted by tidal forces from interactions 
with stars and the disc potential. Therefore P(r) tends 
to zero for light clumps and small distances from the 
Galactic centre. However, more detailed calculations 
show that cuspy centres of the clumps survive the tidal 
disruption. In this case P(r) drops to a constant value 
of .P(O) near the centre of the galaxy [7J. However, P(0) 
is strongly model-dependent, especially it is sensitive 
to the mass distribution of the clumps and the steep- 
ness of the cusp in the centre of the clump. Therefore 



we kept P(0) as a free parameter and left the shape 
of P(r) roughly the same. P(r) is shown in Fig. [2] for 
various values of P(0). The flux can be written as 



@DALclump OC B ■ P(r) ■ p x .clump{l^) dl^ (4) 



where the DM halo profile can be parametrized as 



Px(r)=Po-(-) TTT^ ( 5 ) 





[i + (£) Q l 







In Eq. [5] a is the scale radius of the profile and a, (3 
and 7 are the radial dependencies at r w a, r > a and 
r -C a. The radius ro is the location of the sun and 
po represents the DM density in the solar system, so 
PM = Po- 

The observation of the rotation curves of different dwarf 
galaxies show a flat rotation curve and a constant DM 
density in the Galactic centre, a so-called core. How- 
ever, numerical simulations of galaxy formation prefer 
a Navarro-Frenk- White (NFW) profile [5] which has a 
cusp in the Galactic centre and decreases with I/r 3 
for large radii (a = 1, (3 = 3 and 7=1). Other N- 
body simulations predict that, in addition to these two 
profiles, similar profiles which are even cuspier in the 
centre [9] are possible. 

The question is whether the density profiles /0 XiC / ump (r) 
for the clumpy DM component and p x (r) for the dif- 
fuse DM component are the same or not. According to 
Moore et al. [10] p x is expected to be a cuspy profile, 
while p x . c iump could be a cored profile. 
It was shown in Ref. [5] that the DM annihilation sig- 
nal of a spherical halo profile is not sufficient to re- 
produce the Galactic gamma radiation measured with 
EGRET. To explain the data two ringlike substruc- 
tures in the Galactic disc were needed. A first indi- 
cation of such a substructure was found by Hunter et 
al. in 1997 [11] . The first and inner ring is located at 
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Fig. 3. Longitude distribution of diffuse Galactic gamma radiation above 0.5 GeV with ringlike substructures at 4 kpc 
and f3 kpc 



a radius of 4 kpc. It coincides with a ring of molec- 
ular hydrogen which is located at 4.5 kpc and has a 
width of approximately 2 kpc. The second ring is lo- 
cated at a radius of 13 kpc and most likely resulted 
from the infall of a satellite galaxy. This assumption 
was strengthened by the results of Penarubbia et al. 
[12] who found out that the tidal streams of the dwarf 
galaxy Canis Major form ringlike structures at about 
13 kpc. 

In the z-direction both rings are exponentially dis- 
tributed with a width of 300 pc (830 pc) for the in- 
ner ring (outer ring). However the distribution in r- 
direction is different for the both rings. The inner ring 
has a Gaussian distribution in r-direction with a width 
of 4.3 kpc, while the outer ring has a Gaussian distribu- 
tion with a width of 4.5 kpc to the outside and falling 
of to zero at a radius of 9 kpc by a parametrization of 
an s-shape with two parabolic functions. 
Additional evidence for a ringlike substructure was 



found by Kalberla et al. [13]. They observed the gas 
flaring of molecular hydrogen in the Galactic disc and 
found that they need a ringlike substructure at 13 kpc 
to explain their measured data. This is a very impres- 
sive hint because it is not related to the EGRET ob- 
servations. 

In our analysis a good fit was obtained for a cored 
profile (a = 2, (3 = 4, 7 = 0, a = 5 kpc) for the 
clumpy DM and a NFW profile (a = 1, /3 = 3, 7 = 1, 
a = 100 kpc) for the diffuse DM. The normalization 
parameter po was found to be roughly the same for 
both profiles. In the case of the diffuse DM profile a 
large scale radius was obtained, since the Milky Way is 
located in a local cluster of galaxies which yields large 
amounts of matter at large galactocentric distances. 
The total DM density distributions is the summation 
of all DM components and can be written as 

Px,total Px ' Px: c ^ um P Pinner Ring ~"T~ Pouter Ring (6) 
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Fig. 4. Density of the DM in the Galactic plane. The sum 
of the diffuse and clumpy DM are shown in blue dots, while 
the inner and outer ring are illustrated in green and red 
dots. The position of the sun is marked with a yellow point 



The total mass of the galaxy is constrained by the fact 
that the rotation velocity at the location of the sun is 
v S un = 220 km/ s. The mass of the galaxy within the 
virial radius -R200 is about 10 12 solar masses which is 
in good agreement with Ref . [H] . Note that the diffuse 
component hardly contributes to the DM annihilation 
signal, but significantly to the mass. The density pro- 
files are compared in Fig. O 

The parameters of the halo profiles and the boost fac- 
tor B are varied to minimize the \ 2 function which is 
described in Ref. [2]. The constant boost factor B is 
found to be about 90 and P(0) = 0.7 which agrees well 
with the data and is in the expected range [7]. 
The best fit of the longitude distribution is shown in 
Fig. [3] and the DM density profile with rings is shown 
in Fig. 4. 



3 Summary 

CDM probably consists of WIMPS which are heavy 
and slow particles. If these particles are Majorana par- 
ticles they can annihilate each other and produce Galac- 
tic gamma radiation which can be used to determine 
the density profile of the DM. In this analysis the di- 
rectionality of the DM annihilation flux measured with 
EGRET was used to find a possible DM halo profile. 
After dividing the Galactic gamma ray flux into 4 lati- 
tude and 45 longitude regions the background and the 
DM annihilation signal were fitted to the data for each 
of the 180 bins. The DM annihilation flux is dominated 
by the annihilation flux of the clumpy DM which is 
proportional to p x , c i um p, not p 2 xclump - While the dif- 
fuse DM component has a cuspy NFW profile a shal- 
lower cored distribution was obtained for the clumpy 



component. The DM annihilation flux is dominated 
by the clumpy DM component, but the clumpy com- 
ponent yields a mass below the required mass > 10 12 
solar masses [14] . However, if combined with the dif- 
fuse cuspy NFW profile, both the EGRET data and 
the mass constraint can be fulfilled. In order to take 
the tidal disruption of DM clumps in the vicinity of 
stars into account a survival probability for clumps 
was introduced. Most of the clumps are expected to be 
destroyed near the Galactic centre, although a steep 
cusp may survive. The strong signal observed from 
the Galactic centre yielded a survival probability at 
the centre of P(0) = 0.7. This means that the DM 
clumps are not completely destroyed, which is in good 
agreement with more detailed calculations in Ref. [7|. 
In addition to the DM halo profiles two ringlike sub- 
structure were required at radii of 4 and 13 kpc. The 
halo and ring parameters were obtained by minimizing 
a \ 2 function comparing the flux of the excess from the 
various sky directions with the line-of-sight integral in 
the halo. Figure [3] shows that the halo model fits the 
measured data very well. 

In summary, the EGRET excess of diffuse Galactic 
gamma rays is in good agreement with the expecta- 
tions of a cored clumpy halo component plus a cuspy 
diffuse one. The ringlike substructure, expected from 
the tidal disruption of the nearby Canis Major dwarf 
galaxy, is clearly seen and its heavy mass above 10 10 
solar masses as obtained from the EGRET data, has 
been recently confirmed by the reduced gas flaring in 
this region. 
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